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The  small-amplitude  uniform  magnetization  excitations  in  a  magnetic  nano-pillar  with  two  fer¬ 
romagnetic  layers  coupled  by  the  dipole-dipole  interaction  are  considered.  We  showed  that  the 
dipole  magnetic  field  in  this  system  can  be  fully  described  by  two  diagonal  tensors  for  every  layer 
The  spectrum  of  coupled  oscillations  is  found.  It  consists  of  two  modes,  the  lower  mode  is  quasi- 
antisymmetric  at  small  bias  fields  and  excited  by  current,  the  upper  mode  is  quasi-symmetric  at 
small  bias  fields  and  excited  by  microwave  magnetic  field.  It  is  demonstrated  that  this  spectrum 
can  be  described  by  the  traditional  Kittel  expressions  with  decreased  saturated  magnetization  and 
renormalized  external  magnetic  field.  We  also  demonstrated  that  for  realistic  nano-pillar  parameters 
this  renormalization  does  not  exceed  10  %  for  the  lower  quasi-antisymmetric  mode  at  small  bias 
fields,  and  50  %  for  the  upper  quasi-symmetric  mode  at  small  bias  fields. 

PACS  numbers: 


I.  INTRODUCTION 

The  physics  of  magnetic  multi-layer  systems  with  lay¬ 
ers  of  ferromagnetic  materials  separated  by  non-magnetic 
spacers  with  thickness  in  the  nano-metric  range  is  very 
rich  and  these  systems  have  a  lot  of  practical  applications. 
It  is  well  known  that  the  dipole-dipole  interaction  in  such 
kind  of  systems  plays  crucial  role  due  to  that  the  ferro¬ 
magnetic  layers  are  placed  close  to  each  other  comparing 
to  their  lateral  sizes  and  it  leads  to  substantial  spectrum 
modification  comparing  to  single  layer.  Multi-layer  thing 
films  have  been  existed  and  studied  for  more  than  thirty 
years  while  active  study  of  spin  dynamics  in  multi-layer 
patterned  magnetic  structures  started  relatively  recently 
about  a  decade  ago  when  advances  in  lithography  made 
possible  the  realization  of  patterns  of  nano-metric  mag¬ 
netic  elements  with  controlled  shape  and  dimensions. 

Such  kind  of  systems,  in  which  the  influence  of  dipole 
magnetic  interaction  on  the  spectrum  have  investigated 
recently,  are  ordered  arrays  of  long  magnetic  stripes  with 
micro-metric  width6,7  and  an  array  of  long  magnetic 
stripes  which  are  situated  over  a  continuous  magnetic 
film5 . 

Another  class  of  patterned  system  where  dipole  inter¬ 
action  plays  important  role  are  nano-scale  spin-torque 
oscillators  (STO).  Since  the  first  observations  of  mi¬ 
crowave  generations  by  STO1 1  extensive  experimental  re¬ 
sults  have  been  reported  concerning  the  frequency  mea¬ 
surements  of  this  generation  in  low-amplitude  regime  for 
the  large  variety  of  magnetic  multi-layer  nano-systems. 
In  these  experiments,  the  dependence  of  the  measured 
frequency  on  an  bias  magnetic  field  is  described  well  by 
the  Kittel  expression  for  both  in-plane  and  perpendicu¬ 
larly  magnetized  multi-layers  provided  the  magnetization 
of  the  free  layer  is  essentially  decreased  by  30  to  75  per¬ 
cent  depending  on  a  system4, 11,12,15,18.  The  prevalent 
explanation  of  this  fact  is  the  necessity  of  taking  into 


account  of  a  dipole-dipole  interaction  between  layers. 

Micro-structured  and  nano-structured  ferromagnets 
can  also  be  used  efficiently  as  microwave  absorbers  in 
monolithic  microwave  integrated  circuit  due  to  the  fact 
that  their  FMR  frequency  can  be  widely  tuned  in  the 
range  of  several  tens  of  gigahertz  by  applying  external 
magnetic  fields.  However  constant  application  of  external 
magnetic  fields  is  inconvenient  from  point  of  view  of  low 
electric  power  consuming.  In  order  to  resolve  this  prob¬ 
lem  it  was  proposed  to  control  the  FMR  frequency  by 
using  dipole  magnetic  field  in  patterned  magnetic  multi¬ 
layers16,17. 

In  the  present  paper,  we  theoretically  studied  influence 
of  dipole-dipole  interaction  on  the  spatially-uniform  spin 
wave  modes  in  in-plane  magnetized  two-layered  magnetic 
nano-pillar  shown  on  Fig.  1.  For  simplicity,  we  did 
not  take  into  account  in-plane  anisotropy  (i.e.,  we  as¬ 
sumed  that  the  magnetic  layers  are  patterned  intro  cir¬ 
cular  disks)  and  assumed  that  in  each  disk  the  excited 
mode  has  spatially-uniform  profile.  The  later  assumption 
is  justified  for  small  nano-scale  disks,  in  which  strong  ex¬ 
change  interaction  prohibits  excitation  of  higher  spatial 
modes. 


II.  THE  SYSTEM  OF  EQUATIONS 

The  dynamics  of  coupled  magnetization  oscillations  ex¬ 
ited  by  a  spin  polarized  current  in  a  two  layer  nano-pillar 
is  described  by  two  Landau-Lifshits-Gilbert-Slonczewski 
(LLGS)  equations2, 3,19-21 : 

-gf  -  =7  [H eff,j  X  Mj]  +  Tdj  +  T s,j  +  T (1) 

where  j  =  1,2  and  Mj  is  the  magnetization  of  the  j-th 
layer.  The  first  term  on  the  right-hand  side  of  Eqs.  (1) 
describes  conservative  precession  of  the  magnetization 
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FIG.  1:  (color  online)  Nano-pillar  with  two  ferromagnetic  lay¬ 
ers  coupled  by  the  dipole-dipole  interaction  in  constant  bias 
magnetic  field  H ext,  directed  along  x-axes.  Thicknesses  of  the 
layers  L  1,2  and  the  distance  d  between  them  are  much  smaller 
than  the  layer  radius  R.  The  layers  have  the  same  saturation 
magnetization  M 1  =  M2  =  M  (4-irM  =  8 kOe,  which  is  typical 
for  permalloy). 


vector  Mj  about  the  direction  of  the  effective  magnetic 
field  H eff,j,  where  7  ps  27t  •  2.8GHz/kOe  is  the  modulus 
of  the  gyromagnetic  ratio.  The  structure  of  H ejfj  will 
be  clarified  bellow.  The  second  term 

Toj  =  [Mj  x  [H eff,j  x  Mj]] 

is  dissipative  torque  that  describes  the  energy  dissipation, 
where  a  is  the  damping  constant.  The  third  term 

T.s.j  =  Tn.';  [Mj  x  [Mj  X  Mj/]] 

is  the  Slonczewski-Berger  spin-transfer  torque  that  de¬ 
scribes  interaction  of  the  magnetization  Mj  with  a  spin- 
polarized  current  passing  through  the  layers,  it  depends 
on  orientation  of  the  magnetization  of  the  j  '-layer  ,  where 
j '  =  3  —  j  and  I  is  a  current,  07  =  gefiB /2eMLj  ( g  is 
the  spectroscopic  Lande  factor,  ej  is  a  dimensionless  spin- 
polarization  efficiency,  fis  is  the  Borh  magneton,  e  is  the 
modulus  of  an  electron  charge,  Lj  is  the  thickness  of  the 
j- th  layer).  We  will  choose  positive  direction  of  a  current 
the  way  that  the  sign  ”  is  taken  for  the  first  layer 
and  ”+”  for  the  second  one.  This  means  that  damp¬ 
ing  caused  by  a  spin-polarized  current  will  be  positive 
for  the  first  layer  and  negative  for  the  second  one.  If  we 
neglect  the  magnetization  dynamic  in  the  first  layer  and 
assume  that  its  magnetization  is  fixed  it  will  correspond 
the  well  studied  case  of  magnetization  excitation  only  in 
one  (second)  layer  which  is  called  "free”  in  this  case  (the 
first  layer  is  called  "fixed”).  The  last  term 

T F,j  =  7  [h -(*)  x  Mj] 


is  a  time-dependent  torque  which  describes  action  of  a 
microwave  magnetic  field,  where  h^(f)  =  cos(coextt) 
is  a  microwave  magnetic  field,  Im  is  its  amplitude  and 
ojext  is  its  frequency.  We  will  assume  that  the  vector  lu. 
is  perpendicular  to  static  bias  magnetic  field  Hext- 

Typically,  the  first,  conservative  magnetization  preces¬ 
sion  term  in  the  LLGS  equation  is  much  larger  than  all 
the  other  terms:  dissipation,  spin-transfer  torque,  and 
action  of  a  microwave  magnetic  field  and  we  can  use  the 
following  spin-wave  formalism8,21,  The  solution  of  the 
considered  problem  is  divided  into  two  steps.  At  the  first 
step,  we  consider  the  influence  of  conservative  proces¬ 
sional  torque  and  obtain  spin  wave  modes  by  solving  the 
conservative  Landau-Lifshits  equation  with  precessional 
term  only.  At  the  second  step,  we  analyze  the  influence 
of  damping  and  excitation  torques  on  the  dynamics  of 
previously  found  spin  wave  modes.  It  is  usually  possible 
to  treat  damping  and  excitation  as  small  perturbations, 
which  allows  one  to  obtain  simple  arid  understandable 
analytical  results. 

The  coupled  magnetization  dynamic  which  we  will  con¬ 
sider  can  be  excited  by  either  current  or  microwave  mag¬ 
netic  field. 


III.  THE  DIPOLAR  MAGNETIC  FIELD 
PARAMETERS  FOR  TWO  LAYER 
NANO-PILLAR 


The  Landau-Lifshits  equation  that  describes  non  de¬ 
caying  excitations  in  two  layer  nano-pillar  looks  as  fol¬ 
lows: 

^  =  •>  x  M,1 .  j  =  1, 2.  (2) 

The  effective  magnetic  field  for  j-th  layer,  which 

enters  the  LLGS  equation,  consists  of  the  external  bias 
magnetic  field  H ext  and  magnetodipolar  fields,  created 
by  each  of  the  two  layers: 

2 

He//,j  =  H/./V  +  Hj,fc. 

k= 1 


Here  H77  dipolar  field,  created  by  the  k- th  layer  and  act¬ 
ing  on  the  j- th  one.  The  dipolar  magnetic  field  created  by 
a  layer  depends  on  a  spatial  variable  (it  is  spatially  uni¬ 
form  only  inside  an  uniform  magnetized  ellipse14).  We 
consider  spatially-uniform  spin  wave  excitations  so  the 
field  Hy  /,  should  be  understood  as  the  real  spatially- 
nonuniform  field  H*(r),  which  is  created  by  the  A:-th 
layer,  averaged  over  the  volume  Vj  of  the  j- th  one: 


For  uniformly  magnetized  fc-tli  layer  the  dipolar  field 
H*(r)  can  be  written  as  follows: 


Hfe(r)  =  — 47r.Yfc(r)M*, 


(3) 


3 


where  the  position-dependent  demagnetization  tensor 
Nk(r)  is  given  by  the  expression1  : 

(Awn)  <3*> 

V  J  sP  47r  oXgOXp  J  |  r  —  r | 

vk 

This  expression  for  the  dipolar  field  Hk( r)  is  valid  both 
inside  and  outside  of  the  fc-th  layer.  If  we  calculate 
derivatives  before  taking  the  integral  we  will  get  noth¬ 
ing  else  but  the  ordinary  expression  for  magnetic  field 
produced  by  a  point  dipole.  But  in  this  case,  the  inte¬ 
gral  will  logarithmically  diverge  for  the  diagonal  matrix 
elements  when  the  spatial  variable  r  takes  a  value  inside 
the  k- th  layer.  The  trace  of  the  tensor  Nk{ r)  equals 


Due  to  in-plane  symmetry  of  the  considered  system  of 
two  disks,  demagnetization  tensors  Njk  are  diagonal  in 
the  coordinate  system  represented  on  Fig.  1  and  can  be 
represented  as  follows: 

( Pjk  °  0  \ 

Njk  =  0  Pjk  0  (5) 

v  0  0  Sjk  -  2 pjk  J 

i.e.  each  tensor  depends  only  on  one  parameter  pjk-  In 
the  limit  of  very  thin  disks  pjk  — >  0,  and  dipolar  coupling 
between  layers  vanishes.  By  using  Eqs.  (3a)  and  (3c),  the 
dipolar  magnetic  field  parameters  pjk  can  be  expressed 
in  terms  of  six-fold  multiple  integral  as  follows: 


TrNk(r)  =  -—  /  Ar- - ,=  f  d(r-rk)drk,  (3b) 

47t  J  |r  -  r*  [  J 

14  Vk 

where  (i(r  —  r*)  is  the  Dirac  delta-function. 

By  using  the  above  formulas,  we  can  write  the  dipolar 
fields  Hj-y  as 

H,j,  =  I.T.Vy.M/... 

where 

Njk  =  y  j  •V/,ir  ,),/r;  (3c) 

3  Vj 

are  effective  dimensionless  self-  (for  j  =  k)  and  cross- 
(for  j  ^  k )  demagnetization  tensors,  which  fully  charac¬ 
terize  the  dipole-dipole  interaction  in  our  system.  The 
self-demagnetization  tensor  is  also  called  the  tensor  of 
demagnetizing  coefficients. 

It  should  be  noted  that  the  spatially-nonuniform  self¬ 
demagnetization  tensor  Nj(rj)  for  a  cylinder  (and  cor¬ 
respondingly  for  a  disk)  were  calculated  in  paper10  but 
we  need  spatially-uniform  (averaged  over  a  volume)  both 
self-  and  cross-  demagnetization  tensors. 

Thus,  the  effective  field  for  the  j- th  layer  Heff,j  takes 
a  simple  form: 


"Jt  =  sy  /*■/** 

Vj  Vh 

(6) 

This  integral  can  be  taken  over  an  area  of  the  first  and 
the  second  layers  after  introducing  new  variables  of  inte¬ 
gration,  one  of  which  is  the  difference  of  in-plane  radius- 
vectors  of  the  first  and  the  second  layers  and  the  other 
one  is  in-plane  radius-vector  of  one  of  the  layers.  The 
integral  (6)  can  be  simplified  to 


i  i 


Pjk  ~ 


1  f  .  /  7  r  ( L]  Lk  (1 


0  0 


for  j  ^  k  and 


Pjj  ~ 


(6a) 


h-T  m 


for  j  =  k.  Here 


f(a)  = 


2  +  a 2 


a  \  a 


K - -)  -aE  - 


cr 


H effj  =  Hext  -  4tt  ^  Njk Mfc.  (4) 

fc= l 

The  demagnetization  tensors  Njk  are  symmetrical, 
their  traces  are 

Tr  =  5jk 


where  K  (£)  and  E(£)  are  complete  elliptic  integrals  of 
the  first  and  second  kind. 

The  function  f(a )  has  a  logarithmic  singularity  when 
a  — >  0.  Extracting  this  singularity  and  leaving  terms  up 
to  0(a2  ln(a))  we  have  the  approximate  analytic  expres¬ 
sions  for  the  dipolar  magnetic  field  parameters  in  a  case 
of  thin  disks  (Lip  <C  R  and  d  <C  R): 


where  Sjk  is  the  Kronecker  symbol,  and  and  different 
cross-demagnetization  tensors  are  connected  by  the  ex¬ 
pression: 

VjNjk=VkNkj,  J  / 

These  properties  directly  follow  from  Eqs.  (3a),  (3b), 
(3c)  and  they  are  valid  in  a  general  case,  i.e.  for  a  system 
of  magnetic  bodies  of  arbitrary  shapes. 


„  _  J r_ 

Pjk  ~  2tt  R  '  ° 

- - -  \F(~ 

MLj/R )  L  U 


-  F 


th  f  d  +  Lk 

F(— 


d  +  Lj  \  _ 

R  J 

p  fd  +  LJ  +  Lk 


(6c) 


IV.  THE  SPECTRUM  OF  COUPLED  LINEAR 
EXCITATIONS 
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FIG.  2:  (color  online)  The  dimensionless  dipolar  mag¬ 
netic  field  parameters  pjk,  which  define  self-  and  cross- 
demagnetization  tensors  for  the  system  of  two  thin  disks,  as 
functions  of  the  aspect  ratio  of  the  first  layer  Li/R.  As¬ 
pect  ratio  of  the  second  layer  and  distance  between  layers 
are  fixed:  Li/R,  =  1/10,  d/R  =  1/25.  Solid,  dashed  and 
dotted  lines:  Approximate  analytical  expressions  (6c)  and 
(6d).  Dots:  Exact  expressions  (6a)  and  (6b). 


for  j  ^  k  and 


1  L, 

P»  =  2 nH 


C 


(6d) 


for  j  =  k.  Here  C  =  - 1/2  +  ln(8)  and  F(Q  =  C2  In  Ki¬ 
lt  should  be  noted  that  Eq.  (6d)  can  be  obtained  from 
Eq.  (6c)  by  substituting  Lk  =  Lj  and  d  =  —Lj. 

One  can  see  from  expressions  (6a)  and  (6b)  the  cross- 
demagnetization  tensors  depend  on  aspect  ratio  of  both 
layers  and  relative  distance  between  them:  A/,-/.  = 

Njk(Lj/R,Lk/R,d/R),  the  self-demagnetization  tensors 
depend  on  aspect  ratio  of  correspondent  layer:  Njj  = 
Njj(Lj/R).  When  Lj/R  — >  0  we  have  that  pjj  — >  0 
and  pj'j  — >  0  (j'  =  3  —  j).  In  this  case  the  self¬ 
demagnetization  tensor  of  the  j-th  layer  Njj  passes  into 
the  tensor  of  demagnetizing  coefficients  of  unbounded 
ferromagnetic  film  and  the  cross-demagnetization  tensor 
Nj’j ,  which  describes  the  action  of  the  dipolar  magnetic 
field  created  by  the  j-th  layer  on  the  j'-th  one,  passes 
into  zero  tensor. 

On  Fig  2,  the  dimensionless  dipolar  magnetic  field  pa¬ 
rameters  pjk  are  represented  as  a  function  of  the  as¬ 
pect  ratio  of  the  first  layer  L\/R.  Aspect  ratio  of 
the  second  layer  and  distance  between  layers  are  fixed: 
Lo/R=  1/10,  d/R,  =  1/25. 


We  consider  a  case  of  tangentially  magnetized  nano¬ 
pillar  and  large  enough  constant  bias  magnetic  field, 
when  the  ground  state  of  the  system  is  the  parallel  orien¬ 
tations  of  magnetization  in  both  layers.  We  will  return 
later  to  the  question  on  the  magnitude  of  constant  bias 
magnetic  field  when  the  parallel  orientations  of  magneti¬ 
zation  is  stable.  We  linearize  Eq.  (2)  by  the  substitution: 

M , ( / )  =  M  [x  +  (in,-.  +  c.c.)]  .  (7) 

The  first  term  here  is  the  equilibrium  orientation  of  the 
magnetization  vectors,  whereas  the  second  term  describes 
small-amplitude  (linear)  spin  wave  modes.  The  vectors 
mj  are  orthogonal  to  the  unit  vector  x  in  linear  approx¬ 
imation.  Keeping  only  the  linear  in  m;  terms  yields  the 
following  eigenvalue  problem  for  determination  of  fre¬ 
quencies  u  and  profiles  m j  of  spin  waves: 


iojnij  =7xx 


2 

!  If.  -  1“^ )rn;  • 


2 

+47tM  £  Njk  mfc 

k~i 


(8) 


Condition  of  solvability  of  this  system  determines  fre¬ 
quencies  of  spin  wave  modes: 

(w2  -  w2)  (uj2  -  w2)  =  $intr  (9) 


where 


$int  =  -4pi2P2iw|r  (w2  +  PuP21^m  - 


'2;~,2 


Uh,lUh,2 


4t 4,lW/i,2 


(9a) 


Here  %  =  4wyM,  u>hj  =  7 hj.  The  field  hj  is  the  static 
magnetic  field  acting  on  the  j-th  layer,  i.e.  it  is  external 
field  +  static  dipolar  field  created  by  the  other  layer: 


hj  —  II,  j. j  4tt pjj >  M . 

where  j  '  =  3  —  j.  The  static  dipolar  magnetic  field 
—47 rpjjiM  created  inside  the  j-th  layer  by  the  j'-th  one 
can  reach  a  half  a  kOe  for  realistic  nano-pillar  parameters 
(it  can  be  seen  from  Fig  2).  Frequency  ujj  is  the  preces¬ 
sion  frequency  of  the  j-th  layer  with  taking  into  account 
of  only  static  coupling  to  the  other  layer: 

Zij  =  '•  \Jh j  (h ,  +  4tt[1  -  3pjj]M). 

The  term  <Knt  describes  dynamic  coupling  between  two 
layers  and  depends  on  the  frequency  u>.  When  Lj/R.  —X  0, 
what  means  either  the  thickness  of  the  j-th  layer  tends 
to  zero  or  its  radius  tends  to  infinity,  the  dipolar  mag¬ 
netic  field  produced  by  the  j-th  layer  vanishes  outside  of 
it.  This  lead  to  vanishing  of  dynamic  coupling  between 
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layers  and  function  <I>  ,n(  tends  to  zero.  The  dynamic  cou¬ 
pling  between  layers  is  maximum  when  the  layers  have 
the  same  thickness.  The  term  was  neglected  in  all 
previous  studies. 

The  solutions  of  Eq.  (9)  can  be  approximately  de¬ 
scribed  by  the  traditional  Kittel  expression  with  renor¬ 
malized  external  field  and  saturation  magnetization: 


ah, 2  —  7y  (yHext  —  -Hi, 2)  (^Hext  —  H±*  +  47rMi)2j, 

(10) 

where 

Hi  =  0,  H*  =  47t(pi2  +  p*i)M 

and 


r^Piipii  +  P22P12  +  Zpi2p2i 
P12  +  P21 

_  gPllPl2  +  P'22P2\  ~'2p\2P2\ 

Pl2  +  P21 

The  frequencies  (10)  are  exact  solutions  of  Eq.  (9)  in  two 
limit  cases:  (i)  when  the  layers  have  the  same  thickness 
and  (ii)  when  the  thickness  of  one  of  the  layers  tends  to 
zero. 

The  parameters  Hj  are  found  by  determining  zeros  of 
the  solutions  of  Eq.  (9).  The  expressions  for  renormal¬ 
ized  magnetization  My2  are  the  simplest  ones  that  give 
correct  result  in  both  limit  cases  mentioned  above. 

For  symmetrical  system  (both  layers  have  identical 
thicknesses),  mode  with  the  frequency  Ui  corresponds 
to  symmetric  excitations  (the  magnetization  precession 
in  both  layers  have  the  same  phases),  whereas  mode 
with  the  frequency  cu2  describes  antisymmetric  excita¬ 
tions  (opposite  phases  of  precession  in  both  layers).  For 
different  thicknesses  of  the  layers  one  can  still  classify 
the  mode  Ui  as  quasi-symmetric,  and  inode  w2  as  quasi- 
antisymmetric. 

The  frequencies  of  coupled  linear  excitations  are  rep¬ 
resented  on  Fig.  3.  One  can  see  that  expressions  (10) 
fit  the  solution  of  secular  Eq.  (9)  with  high  accuracy  at 
small  bias  magnetic  fields.  When 

Hext  =  H2(4ttM2  -  Ho)/( 4n(Mi  -  M2)  -  2 Ho) 

the  curves  that  correspond  to  approximate  expressions 
(10)  intersect  while  the  curves  that  correspond  to  the 
solution  of  secular  Eq.  (9)  do  not.  The  modes  described 
by  secular  Eq.  (9)  change  their  symmetry  near  this  point, 
while  the  modes  described  by  approximate  expressions 
(10)  retain  their  symmetry.  On  Fig  4,  the  maximum 
and  minimum  angle  between  the  vectors  mi  and  m2, 
which  describe  the  coupled  linear  spin  wave  excitation  in 
the  first  and  the  second  layer  respectively,  for  frequencies 
that  are  the  solutions  of  Eq.  (9)  are  represented.  While 
calculating  we  put:  L\/R  =  1/5 ,Lo/R  =  1/10, d/R.  = 
1/25. 


FIG.  3:  (color  online)  (a)  The  frequencies  of  linear  cou¬ 
pled  excitations  in  a  two  layer  nano-pillar  as  a  function 
of  bias  magnetic  field.  Solid  lines:  The  exact  solutions 
of  secular  Eq.  (9).  Triangles  and  dots:  The  fre¬ 
quencies  ici/27t  and  0:2/277  respectively  (approximate  ex¬ 
pressions  (10)).  Dashed  and  dotted  lines:  The  tradi¬ 
tional  Kittel  expressions  for  in-plane  magnetized  ferromag¬ 
netic  film13:  w/v.i  =  ^i/Hext{Hext  +47 rtf)  and  coa',2  = 

7 ij {He.ri  —  Ho)(Hext  -  H-2  +  47 tM)  respectively.  (b)  The 
frequencies  (a;  —  yHext) /27r  as  a  function  of  bias  magnetic 
field.  Solid  lines:  The  exact  solutions  of  secular  Eq.  (9). 
Triangles  and  dots:  (on  — 'jHtxt) /2tt,  (cj2  ~  jHext) /2tt  re¬ 
spectively  (approximate  expressions  (10)).  While  calculating 
we  put:  Li/R  =  1/5 ,L-2/R  =  1/10, d/R  =  1/25. 


The  parallel  magnetization  ground  state  is  stable  when 
external  bias  magnetic  field  is  larger  than  H2  (Ho  = 
0.7 kOe  for  the  spectrum  represented  on  Fig.  3)  otherwise 
one  of  the  frequencies  is  complex  what  leads  to  instability 
of  the  ground  state. 

As  one  can  see  from  Fig.  3  the  both  modes  lie  lower 
corresponding  Kittel  modes,  that  means  the  both  mode 
correspond  to  effective  decrease  of  magnetization  and  the 
decrease  for  quasi-symmetric  mode  is  larger.  The  rela- 
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L-z/R  =  1/10,  d/R  =  1/25.  One  can  see  from  Fig.  5 
when  the  difference  in  layer  thickness  increases  (it  means 
the  coupling  between  layers  get  weaker)  the  effective  de¬ 
crease  of  magnetization  occurs  mainly  due  to  correspond¬ 
ing  self-demagnetization  tensors  (the  tensors  of  demagne¬ 
tizing  coefficients).  As  one  can  see  from  Fig.  5  the  effec¬ 
tive  decrease  of  magnetization  for  upper  quasi-symmetric 
mode  at  small  bias  fields  can  reach  50  %  for  still  realistic 
nano-pillar  parameters  (due  to  self-demagnetization  ten¬ 
sor  mainly  )  and  effective  decrease  of  magnetization  for 
lower  quasi-antisymmetric  mode  can  only  be  about  10%. 

V.  THE  EXCITATION  OF  COUPLED  LINEAR 
SPIN  WAVE  MODES  BY  CURRENT  AND 
MICROWAVE  MAGNETIC  FIELD 


FIG.  4:  (color  online)  The  maximum  and  minimum  angle 
between  the  vectors  mi  and  m2,  which  describe  the  coupled 
linear  spin  wave  excitation  in  the  first  and  the  second  layer 
respectively,  for  frequencies  that  are  the  solutions  of  Eq.  (9) 
as  a  function  of  a  bias  magnetic  field.  Solid  red  lines:  for 
lower  mode  on  Fig  3.  Dashed  blue  lines:  for  upper  mode 
on  Fig  3.  While  calculating  we  put:  Li/R  =  l/b,L2/R  = 
1/10,  d/R  =  1/25. 


Aspect  ratio  of  the  first  layer  L\  /  R 


FIG.  5:  (color  online)  Solid  lines:  Relative  decrease  of  effec¬ 
tive  magnetization  (Mi, 2  —  M)/M  for  both  quasi-symmetric 
and  quasi-antisymmetric  modes  with  respect  to  the  mag¬ 
netization  of  unbounded  ferromagnetic  film.  Dotted  and 
dashed  lines:  The  effective  decrease  of  magnetization  for  the 
first  and  the  second  isolated  (uncoupled)  layers,  respectively. 
While  calculating  we  fixed  the  thickness  of  the  second  layer 
and  the  distance  between  layers:  Li/R  =  1/10,  d/R  =  1/25. 


five  decrease  of  magnetization  (Mi, 2  —  M)/M  with  re¬ 
spect  to  the  magnetization  of  unbounded  ferromagnetic 
film  as  a  function  of  aspect  ratio  of  the  first  layer  is  rep¬ 
resented  on  Fig.  5.  While  calculating  we  fixed  the  thick¬ 
ness  of  the  second  layer  and  the  distance  between  layers: 


The  eigenvalue  problem  (8)  give  us  four  different  fre¬ 
quencies  (two  positive  and  two  negative).  It  means  all 
eigenvectors  {m^',  m^'}  (index  v  =  1,4  numerate  dif¬ 
ferent  eigenvalues),  which  describe  the  coupled  linear 
spin  wave  excitation  in  the  first  and  second  layers,  are  lin¬ 
ear  independent  and  they  form  the  linear  space  C2  x  C2 . 
They  are  orthogonal  with  respect  to  the  inner  product9: 

2 

-i  ^2  Vjx  ■  [mb')*  x  m19'  ;]  =  Nv5vvi, 

3= 1 

where  N„  is  a  norm  (defined  by  this  expression)  and  <W 
is  the  Kronecker  symbol. 

We  can  expand  the  solution  of  linearized  Eq.  (1)  in 
a  full  set  of  functions  defined  by  eigenvalue  problem  (8). 
After  substitution 


M  j(t)=M 


4 

x  +  y^(m (pcv{t)  +  c.c.) 

V=1 


(ID 


into  Eq.  (1)  linearizing  it  and  taking  the  inner  product 
with  we  get  the  system  of  equations  that  is 

equivalent  to  linearized  Eq.  (1): 


^  ^  —  ^2  i^vdpp'  +  i[otQi/v'  —  IPVv'])cp'  + 

''  p'=i 

T  S„  co s ( ce,..,, /  f ) ,  (12) 

where  is  a  frequency  of  linear  coupled  excitation  de¬ 
fined  by  a  secular  Eq.  (9)  and 


Qw<  =-jf52  Vi  ■  <  ’)  (12a) 

"  3= 1 

are  matrix  elements  which  define  damping  caused  by  the 
damping  torque, 

PVp'  =-^~J2  G3Vi  x  ‘  [“i^*  x  (m2  ]  ~  *)]  (12b) 
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are  matrix  elements  which  define  negative  damping 
caused  by  the  Slonczwski-Berger  torque.  Both  matrixes 
Q  arid  P  are  real.  And 

^  =  (12c) 
v  i— 1 

we  assumed  that  the  vector  of  amplitude  of  microwave 
magnetic  field  h^,  is  perpendicular  to  constant  bias  mag¬ 
netic  field  Hea.(. 

We  will  use  the  fact  that  conservative  magnetization 
precession  term  in  the  LLGS  equation  is  much  larger 
than  dissipation,  spin-transfer  torques,  and  action  of  a 
microwave  magnetic  field  and  we  will  treat  them  as  small 
perturbation  so  we  will  solve  Eq.  (12)  using  perturbation 
theory. 

We  will  consider  the  action  of  current  and  microwave 
magnetic  field  separately.  We  will  consider  the  excitation 
of  coupled  linear  spin  wave  modes  by  current  first  so  we 
consider  Eq.  (12)  with  S„  =  0.  We  will  use  the  following 
substitution  c„  =  elQvt .  We  have 

4 

F  ([f^  -  u)v]5vv>  -  i[aQvv>  -  JP„„/])<v  =  0,  (13) 

v'=l 

here  fi„  is  a  complex  frequency,  the  real  part  of  which  is 
the  frequency  of  coupled  linear  excitation  and  imaginary 
part  is  the  damping.  The  imaginary  part  depends  on  a 
current  I  and  the  Gilbert  constant  a.  By  equating  it 
to  zero  we  get  the  threshold  current  Ith  for  excitation  of 
coupled  linear  spin  wave  modes. 

The  first  order  of  the  perturbation  theory  give  us: 

fi;,  =  +  i[aQvv  -  IPVV ].  (14) 

One  can  see  the  damping  is  linear  with  respect  to  cur¬ 
rent  I  and  the  Gilbert  constant  a  in  the  first  order  of 
the  perturbation  theory.  The  second  order  of  the  pertur¬ 
bation  theory  is  real  and  give  us  the  adjustment  to  the 
frequency  of  the  coupled  linear  excitation. 

The  threshold  current  for  excitation  of  coupled  linear 
spin  wave  modes  as  a  function  of  a  bias  magnetic  field  is 
represented  on  Fig.  6.  One  can  see  that  current  excites 
lower  (quasi-antisymmetric)  mode  at  small  bias  fields 
and  upper  mode  (which  passes  from  quasi-symmetric  to 
quasi-antisymmetric  with  bias  field  increase)  at  large  bias 
fields.  It  should  be  noted  that  both  frequencies  (see 
Fig.  3)  are  close  to  each  other,  do  not  intersect  at  large 
bias  fields  and  so  it  is  not  very  important  practically 
which  mode  is  excited  at  these  bias  fields. 

By  using  expression  (14),  the  expression  for  threshold 
current  at  small  bias  fields  can  be  written  as  follows: 

hh  =  a— ,  (lo) 

4  22 

where  Q 22  and  P22  are  the  matrix  elements  (12a)  and 
(12b)  calculated  on  the  eigenfunctions  which  correspond 


FIG.  6:  (color  online)  The  threshold  current  for  excitation 
of  coupled  linear  spin  wave  modes  represented  on  Fig.  3  as  a 
function  of  a  bias  magnetic  field.  While  calculating  we  put: 
Li/R  =  1/5,  L2/R  =  1/10,  d/R  =  1/25,  R  =  50 nm,a  = 
0.01,  si  =  £2  =  1.  Solid  and  dashed  lines:  The  thresh¬ 
old  current  for  excitation  of  lower  mode  (quasi-antisymmetric 
one  at  small  bias  external  fields)  and  upper  mode  (quasi- 
symmetric  one  at  small  bias  external  fields),  respectively,  ob¬ 
tained  from  approximate  Eq.  (14).  Dots  and  triangles: 
The  same  threshold  current  obtained  from  Eq.  (13). 


Aspect  ratio  of  the  second  layer  L2/R 


FIG.  7:  (color  online)  The  threshold  current  for  excitation  of 
lower  mode,  which  is  quasi-antisymmetric  one  at  small  bias 
fields,  as  a  function  of  aspect  ratio  of  the  second  layer  Li/R 
when  external  bias  magnetic  field  Hext  =  3 kOe.  While  cal¬ 
culating  we  put:  Li/R  =  1/10, d/R  =  1/25,1?  =  50 urn, a  = 
0.01,  £1  =£2  =  1.  Solid  line:  approximate  expression  (15). 
Dashed  line:  approximate  expressions  (15a)  and  (15b). 
Dots:  the  threshold  current  obtained  from  Eq.  (13). 


to  lower  mode,  which  is  quasi-antisymmetric  at  small  bias 
fields. 

The  matrix  elements  Qvv'  define  damping  and  depend 
on  layer  sizes  weakly.  Matrix  elements  Pvv'  =  0  for  a 
symmetrical  nano-pillar  (when  the  layers  have  the  same 
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thickness).  It  means  current  can  excite  neither  mode  in 
this  case. 


Bias  magnetic  field  Hext  (kOe) 

FIG.  8:  (color  online)  The  coefficients  | S\,  which  show  the 
amplitude  of  excitation  of  coupled  spin  wave  modes  by  mi¬ 
crowave  magnetic  field,  as  a  function  of  a  bias  magnetic  field. 
While  calculating  we  put:  Li/R  =  1/10,  d/R  =  1/25  and 
h.,  =  {0. 0.  IkOe.  ()}.  Solid  and  dashed  lines:  The  coeffi¬ 
cients  151  for  lower  mode  (quasi-antisymmetric  one  at  small 
bias  external  fields)  and  upper  mode  (quasi-symmetric  one  at 
small  bias  external  fields),  respectively,  obtained  from  expres¬ 
sion  (12c)  .  Dotted  and  dashed  dotted  lines:  The  same 
coefficients  |5|  obtained  from  approximate  expressions  (16a) 
and  (16b). 

We  can  decompose  the  matrix  elements  Q 22  and  P22 
with  respect  to  the  ratio  L2/Li  and  restrict  ourselves  by 
zero  and  first  nonzero  approximations  respectively.  We 
have 

Q22  =  j(H  +  —  2>pn  —  P21)),  (15a) 

47rM/>2i(— 3pn  +  A3A2)  —  H(—3pn  +  AiA2) 

100  — - X 

2P21A2(4ttMA4- H) 

x  (i  -  ,  (15b) 

where  Ai,  A2,  A3,  and  A4  are  dimensionless  functions: 

Ai  =  3pn  —p2i,  A2  =  C-HLi/R),  A3  =  2  — 3pn  —  p2\, 

A4  =  1  —  3pn  4-  P21  • 

The  threshold  current  for  excitation  of  lower  mode, 
which  is  quasi-antisymmetric  one  at  small  bias  fields,  as 
a  function  of  aspect  ratio  of  the  second  layer  L2 /R  when 
external  bias  magnetic  field  Hext  =  3 kOe  is  represented 
on  Fig.  7. 

We  will  consider  the  excitation  of  coupled  linear  spin 
wave  modes  by  microwave  magnetic  field  so  we  consider 
the  Eq.  (12)  with  I  =  0.  The  mode  for  which  the  abso¬ 
lute  value  of  the  coefficient  S  is  larger  has  larger  ampli¬ 
tude  and  it  will  be  excited  by  microwave  magnetic  field. 

The  coefficients  151  for  both  modes  are  represented  on 
Fig.  8.  One  can  see  that  microwave  magnetic  field  excites 


upper  (quasi-symmetric)  mode  at  small  bias  fields  and 
lower  mode  (which  passes  from  quasi-antisymmetric  to 
quasi-symmetric  with  bias  field  increase)  at  large  bias 
fields. 

The  same  way  as  we  get  approximate  expressions  for 
Q22  and  P22  we  can  obtain  approximate  expressions  for 
Si  and  S2  calculated  on  the  eigenfunctions  which  cor¬ 
respond  to  the  upper  and  lower  mode  respectively.  We 
have 

Si  =  ~  \h~ 7+ 

3y  h^(H  —  8irM  p2i)(pi4  —  A5A2)  /  L2  \ 

8p2  \A2(H  —  4n M  A.) )  \  Lj  ’ 

(16a) 


_  3y h^(H  -  87rAfp2i)(pn  -  A5A2)  /  _  LA 
-  8p21A2(H -4nMA4)  V  A )  ' 

(16b) 

where  A5  =  pu  —  p2\.  We  assumed  the  vector  is 
directed  along  y- axis. 


VI.  CONCLUSIONS 

We  showed  the  dipole  magnetic  field  in  a  two  layer 
nano-pillar  can  be  fully  described  by  two  diagonal  ten¬ 
sors  for  every  layer:  one  of  them  is  the  usual  tensor  of 
demagnetizing  coefficients  (self-demagnetization  tensor), 
the  other  one  is  the  cross-demagnetization  tensor  which 
describes  the  dipole-dipole  interaction  between  different 
layers.  Every  tensor  is  fully  characterized  by  only  one 
dimensionless  parameter  which  depends  only  on  relative 
geometrical  sizes  of  the  nano-pillar  system. 

The  spectrum  of  coupled  oscillations  of  a  two-layer 
nano-pillar  consists  of  two  modes  with  different  fre¬ 
quencies.  We  showed  that  the  lower  mode  is  quasi- 
antisymmetric  at  small  bias  fields  and  excited  by  cur¬ 
rent,  the  upper  mode  is  quasi-symmetric  at  small  bias 
fields  and  excited  by  microwave  magnetic  field. 

We  showed  that  the  frequencies  of  coupled  spin  wave 
modes  can  be  described  by  the  traditional  Kittel  expres¬ 
sion  with  renormalized  bias  field  and  reduced  saturation 
magnetization.  The  dipole  magnetic  field  leads  to  effec¬ 
tive  decrease  of  magnetization  for  both  modes  but  this 
decrease  is  different  for  different  modes.  It  does  not  ex¬ 
ceed  10  %  for  the  lower  quasi-antisymmetric  mode  at 
small  bias  fields  (which  is  excited  by  spin-polarized  cur¬ 
rent),  and  can  be  50  %  for  the  upper  quasi-symmetric 
mode  at  small  bias  fields  (which  is  excited  by  microwave 
magnetic  field)  For  still  realistic  nano-pillar  parameters. 
But  it  should  be  noted  these  50  %  occur  mainly  due  to 
the  tensor  of  demagnetizing  coefficients  (not  due  to  the 
dynamic  dipole-dipole  interaction  between  layers). 
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